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Abstract 

We construct classes of stochastic differential equations with fluc- 
tuating friction forces that generate a dynamics correctly described 
by Tsallis statistics and nonextensive statistical mechanics. These 
systems generalize the way in which ordinary Langevin equations un- 
derly ordinary statistical mechanics to the more general nonextensive 
case. As a main example, we construct a dynamical model of velocity 
fluctuations in a turbulent flow, which generates probability densities 
that very well fit experimentally measured probability densities in Eu- 
lerian and Lagrangian turbulence. Our approach provides a dynamical 
reason why many physical systems with fluctuations in temperature 
or energy dissipation rate are correctly described by Tsallis statistics. 
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Recently there has been considerable interest in the formalism of nonex- 
tensive statistical mechanics (NESM) as introduced by Tsallis [|l| and further 
developed by many others (e.g. @-@]). In the mean time there is growing 
evidence that the formalism, rather than being just a theoretical construc- 
tion, is of relevance to many complex physical systems. Applications in 
various areas have been reported, mainly for systems with either long-range 
interactions 0-0, multifractal behaviour |§ |J, or fluctuations of temper- 
ature or energy dissipation rate ||10||— 1|14||. A recent interesting application 



of the formalism is that to fully developed turbulence || [TT], [T2[]. Precision 
measurements of probability density functions (pdfs) of longitudinal veloc- 
ity differences in high-Reynolds number turbulent Couette- Taylor flows are 
found to agree quite perfectly with analytic formulas of pdfs as predicted by 
NESM pj. 

Despite this apparent success of the nonextensive approach, still the ques- 
tion remains why in many cases (such as the above turbulent flow) NESM 
works so well. To answer this question, let us first go back to ordinary sta- 
tistical mechanics and just consider a very simple well known example, the 



Brownian particle fl5 |. Its velocity u satisfies the linear Langevin equation 

u = — 7M + crL(t), (1) 

where L(t) is Gaussian white noise, 7 > is a friction constant, and a 
describes the strength of the noise. The stationary probability density of u 
is Gaussian with average and variance j3~ l , where where (3 = \ can be 
identified with the inverse temperature of ordinary statistical mechanics (we 
assume that the Brownian particle has mass 1). 

The above simple situation completely changes if one allows the param- 
eters 7 and a in the stochastic differential equation (SDE) to fluctuate as 
well. To be specific, let us assume that either 7 or o or both fluctuate in 
such a way that (3 = 7/c 2 is x 2 -distributed with degree n. This means the 
probability density of (3 is given by 

n 

w) = FT!)fe}>M-if} (2) 

The x 2 distribution is a typical distribution that naturally arises in many 
circumstances. For example, consider n independent Gaussian random vari- 
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ables Xi, i = 1, . . . , n with average 0. If (3 is given by the sum 

n 

p--=Y, x " ( 3 ) 

i=l 

then it has the pdf (0). The average is given by 

POO 

((3) = n(X 2 ) = / Pf((3)dp = (3 (4) 
Jo 

and the variance by 

(P 2 ) ~ P 2 o = hi (5) 

(see also pf). 

Now assume that the time scale on which (3 fluctuates is much larger 
than the typical time scale of order r y~ 1 that the Langevin system ([I]) needs 
to reach equilibrium. In this case one obtains for the conditional probability 
p{u\f3) (i.e. the probability of u given some value of (3) 

^ . f 1 , 2 



pW) = V^ eXP l 2^ J ' (6) 

for the joint probability p(u, j3) (i.e. the probability to observe both a certain 
value of u and a certain value of j3) 

P ( Uj f3)=p(u\P)f(P) (7) 

and for the marginal probability p(u) (i.e. the probability to observe a certain 
value of u no matter what f3 is) 



p(u) = j p{u\(3)f{(3)d(3. (8) 
The integral (|8]) is easily evaluated and one obtains 

Hence the SDE (|l]) with x 2 -distributed (3 = 7/cx 2 generates the generalized 
canonical distributions of NESM |J 

p{u) ^ (10) 

q-l 



(l + f /% - l)u< 



provided the following identifications are made. 



1 n 1 , 2 tiis 
77 + 77 q = l + (11) 



q-1 2 2 n+1 

kq-l)P = - ^ = T^-/3o- (12) 

2 n 3 — g 

We already see from this simple example that the physical inverse tempera- 
ture (3 = (P) not necessarily coincides with the inverse temperature param- 
eter (3 used in the nonextensive formalism (see also H for related results). 

More generally, we may also consider nonlinear Langevin equations of the 
form 

u = -jF(u) + aL(t) (13) 

where F(u) = —J^V(u) is a nonlinear forcing. To be specific, let us assume 
that V(u) = C\u\ 2a is a power-law potential. The SDE (JT3|) then generates 
the conditional pdf 

p(u\P) = — ^ (C0)& exp {-PC\u\ 2a } (14) 

and for the marginal distributions p(u) = f p(u\/3)f((3)d/3 we obtain after a 
short calculation 

P{u) = ^ \ -, (15) 

where 

r W 

Z" 1 = a f C( ff - 1)^) *" • — V V , N (16) 



\2a) \q-l 2a J 



and 



an + 1 
2<y 

P = y— 9 A)- (18) 

1 + 2a — q 

To generalize to iV particles in d space dimensions, we may consider 
coupled systems of SDEs with fluctuating friction forces, as given by 

Ui = -jiFi («!,... ,u N ) +<T i L i (t) i = l,...,N (19) 
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Suppose that a potential V(u[, . . . , tT N ) exists for this problem such that 
Fi = -^-V. Moreover, assume that all fii = ^ fluctuate in the same way, 

i.e. are given by the same fluctuating \ 2 distributed random variable fii = fi. 
One then has for the conditional probability 

p(u h ...,u N \fi) = — exp {-fiV{u u . . . , u N )}, (20) 

where Z(fi) = J du\ . . . du^ e~^ v is the partition function of ordinary statis- 
tical mechanics. Suppose that Z(fi) is of the form 

Z(fi) ~ p x e~P y , (21) 

then integration over the fluctuating fi leads to marginal distributions of the 
form ^ 

p(tTi, ...,vT N ) — , (22) 

(1 + ^-^1,...,%))^ 

i.e. the generalized canonical distributions of NESM with 

« = 1 + ^ < 23) 

and 

n = ^° (oA\ 

Eq. (|22| ) is correct if all particles see the same fluctuating fi at the same 
time — an assumption that can only be true for a very dense and concentrated 
system of particles. In many physical applications, the various particles will 
be dilute and only weakly interacting. Hence in this case fi is expected to 
fluctuate spatially in such a way that the local inverse temperature fii sur- 
rounding one particle i is almost independent from the local fij surrounding 
another particle j. Moreover, the potential is approximately just a sum of 
single-particle potentials V(wi, • • • , Wjv) = YuiLi K(wj). In this case integra- 
tion over all fii leads to marginal densities of the form 

N 1 

p(ux,...,u N ) ~TT — , (25) 

l\ (1 + fi(q - l)V s (u t ))- 

i.e. the iV-particle nonextensive system reduces to products of 1-particle 
nonextensive systems (this type of factorization was e.g. successfully used 
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in The truth of what the correct nonextensive thermodynamic descrip- 

tion is will often lie inbetween the two extreme cases (|22|) and (p5|). 

Let us now come to our main physical example, namely fully developed 
turbulence. Let u in eq. (p~3|) represent a local velocity difference in a fully 
developed turbulent flow as measured on a certain scale r. We define 

P = er (26) 

where e is the (fluctuating) energy dissipation rate (averaged over r 3 ) and 
r is a typical time scale during which energy is transferred. Both e and 
t can fluctuate, and we assume that er is x 2_ distributed. For power-law 
friction forces the SDE ( |T3| ) generates the stationary pdf flip]). In Fig. 1 



this theoretical distribution is compared with experimental measurements 
in two turbulence experiments, performed on two very different scales. All 
distributions have been rescaled to variance 1. Apparently, there is very good 
coincidence between experimental and theoretical curves, thus indicating that 
our simple model assumptions are a good approximation of the true turbulent 
statistics. 

Generally, in turbulent systems the entropic index q is observed to de- 



crease with increasing r (see |L2 for precision measurements of q(r)). At the 



smallest scale we expect from the definition of the energy dissipation rate 

(27) 




[y: kinematic viscosity, v. velocity) that the independent Gaussian random 
variables Xi in eq. @ are given by 

dv 



Xi = v 7 ^^-;, (28) 

and that there are indeed 3 of them, due to the 3 spatial dimensions. This 
means n = 3 or, using eq. (0), q ~ § if a ~ 1. This is indeed confirmed by 
the fit of the small-scale data of the Bodenschatz group in Fig. 1, yielding 
q = 1.49 and an a as given by eq. (|17|). 
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Fig. 1 Histogram of longitudinal velocity differences as measured by 



Swinney et al. [|12|, in a turbulent Couette Taylor flow with Reynolds 
number R\ = 262 at scale r = 116^ (solid line), where 7] is the Kolmogorov 
length. The experimental data are very well fitted by the analytic formula 
(|j~5| ) with q = 1.10 and a = 0.90 (dashed line). The square data points are 
a histogram of the acceleration (= velocity difference on a very small time 
scale) of a Lagrangian test particle as measured by Bodenschatz et al. for 
R x = 200 0. These data are well fitted by (0) with q = 1.49 and a = 0.92 
(dotted line). 
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